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ABSTRACT 



Aims. We have obtained and analyzed Chandra ACIS-S observations of the strongly accreting classical T Tauri star DG Tau. Our 
principal goals are to map the immediate environment of the star to characterize possible extended X-rays formed in the jet, and to 
re- visit the anomalous, doubly absorbed X-ray spectrum of DG Tau itself. 

Methods. We combine our new ACIS-S data with a data set previously obtained. The data are superimposed to obtain flux and 

hardness images. Separate X-ray spectra are extracted for DG Tau and areas outside its point spread fu nction. 

Results. We detect a prominent X-ray jet at a position angle of PA w 225 deg (tentatively suggested by Giidel et al. 2005), coincident 
with the optical jet axis. We also identify a counter jet at PA = 45 deg. The X-ray jets are detected out to a distance of = 5" from the 
star, their sources being extended at the ACIS-S resolution. The jet spectra are soft, with a best-fit electron temperature of 3.4 MK. We 
find evidence for excess absorption of the counter jet. The spectrum of the DG Tau point source shows two components with largely 
different temperatures and absorption column densities. 

Conclusions. The similar temperatures and small absorbing gas columns of the jet sources and the soft component of the "stellar" 
source suggest that these sources are related, produced either by shocks or by magnetic heating in the jets. Cooling estimates suggest 
that the pressure in the hot gas contributes to jet expansion. The hard "stellar" component, on the other hand, is associated with a 
stellar corona or magnetosphere. The excessive photoelectric absorption of this component suggests the presence of dust-depleted 
accretion streams above coronal magnetic fields. 
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1. Introduction 

Pre-main sequence stars show various signs of accretion and out- 
flow, such as stella r winds inducing high mass-loss rates (e.g., 
iDupree et all l2005t Uohns-Krull & Herczegl I20071; iKwan et alj 
2007) , molecul ar outflows observed in molecular lines (e.g., 
Bachiller 1996), and accompa nying optical (e.g., Hirth et al. 
19971; lEisloffel & Mundtl 1 1998b and radio jets (e.g., [Anglada 
1995b . The most evident manifestation of outflows are the op- 



tically visible jets and their associated Herbig-Haro (HH) ob- 
jects at distances up to several arcminutes from the star. These 
structures are excited by internal shocks or in regions where 
the fast mass stream encounters the interstellar medium and 
shock-ionizes the gas ( for a review of Herbig-Haro flows, see 
Reipurth & Bally 2001). Under ideal circumstances (low extinc- 
tion, strong ionization), o ptical jets can be ident ified at distances 
as close as 07 1 to the star (Bacciot ti et afl,l2002l) . The same com- 
pact jets are also routinely detected at radio wavelengths, where 
the emission mechani sm is thought to be bremsstrahlun g from 
the shock-heated gas dRodri guezl l 1 9951; IXngladal 1 1 9951) . Radio 
brightness temperatures suggest overall gas temperatures of or- 
der 10 4 K. This picture is ambiguous, however, as a number 
of non-thermal jets have been suggested from radio polariza- 



tion or synchrot ron-like spectral shapes (e.g. JCuriel et aljfl 993h 
iRav et ai]|1997l) . Magnetic fields may thus play a role not only 
in launching the jets, but in their propagation as well. 

Outflow processes are prone to producing X-rays, given that 
shocks with shock jump velocities of order several hundred 
km s" 1 are possible. The relevant theo ry and a simple model have 
been discussed bv lRaga et al.1 d2002l) . The strong-shock temper- 
ature can be expressed as T * 1.5 x 10 5 Vj ()0 K (for fully ionized 
gas) where vioo is the shock speed relative to a target in units of 
100 km s -1 . Jet speeds are typically of order v = 300-500 km s _1 
(lEisloffel & Mundl 1 1998tlAngladal.il 9951 IBallv et al.Ll2003h . in 
principle allowing for shock speeds of similar magnitude. If a 
flow shocks a standing medium at 400 km s , then T » 2.4 MK. 



Faint, soft X-ray emission has been detected 
a few protostellar HH objects (Pravdo e t all 12001 



iPravdo & Tsubol 12001 iFavata et all Hool iBallv et all 



from 
2001: 



2003] 

Tsuiimoto et all 120041 iGrosso et all 120061) . IBallv et all (2003) 
used a Chandra observation to show that X-rays form within 
an arcsecond of the protostar L1551 IRS-5 while the star itself 
is too heavily obscured to be detected. As this example illus- 
trates, the jet-launching region of powerful protostellar jets is 
often inaccessible to optical, near-infrared or X-ray studies due 
to excessive absorption. However, a class of strongly accreting, 
optically revealed classical T Tauri stars (CTTS) also exhibits so- 
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Table 1. Observations 



Previous observations: 



Instrument 


Chandra ACIS-S 


XMM-Newton EPIC 




ObsID 


4487 


0203540201 




Start time (UT) 


2004-01-11 02:58:51 


2004-08-17 06:08:10 




End time (UT) 


2004-01-11 11:52:21 


2004-08-17 17:32:46 




Exposure time (s) 


29717 


41076 




New observations: 


Instrument 


Chandra ACIS-S 


Chandra ACIS-S 


Chandra ACIS-S 


ObsID 


6409 


7247 


7246 


Start time (UT) 


2005-12-15 11:03:14 


2005-12-17 09:37:05 


2006-04-12 17:13:23 


End time (UT) 


2005-12-15 16:17:01 


2005-12-17 14:34:37 


2006-04-13 01:37:19 


Exposure time (s) 


16252 


15946 


27811 



called micro-jets visible in optical lines dHirth et al.L [l997). with 
flow speeds similar to protostellar jets. CTTS micro-jets have the 
unique advantage that they can - in principle - be followed down 
to the acceleration region close to the sta r both in the o ptical and 
in X-rays. For example, iBacciotti et al.l d2000l [2002) used the 
HST to trace the jet of the CTTS DG Tau to within 071 of the 
star. 

DG Tau is a most outstanding T Tauri star for X-ray studies. 
A Chandra high-resolution X-ray image has shown tentative evi- 
dence for the p resence of fa i nt X-ra ys along the optic al jet. Both 
XMM -Newton dGiidel et all l2007bl) and Chandra dGiidel et all 
120051) low-resolution CCD spectra of DG Tau are anomalous, 
showing a "two-absorber X-ray" (TAX) spectrum in which two 
independent X-ray components are each subject to different ab- 
sorption column densities. 

To study DG Tau's X-ray emission further, we obtained new 
Chandra observations that, together with the previous observa- 
tions, result in three times more Chandra exposure than analyzed 
before. The present paper describes the new observations, puts 
them into a context of the previous results, and discusses some 
tentative models for the jet X-ray emission. 



11", with velocities of several 100 km s 1 (e.g. , Lavallev et al.l 
119971: lEisloffel&Mundtlll998t iDougados et afl |2000l). Optical 
evidence for a counte r jet has been reported dLavallev et all 
fT997tlPvo etall 120031 

2.2. Previous X-ray observations 

We have previously obtained two short X-ray observations of 
DG Tau (Table ED- Our previous Chandra ACIS-S observation 
dGiidel etall [2005) showed tentative indications of very faint, 
soft emission along the optical forward jet out to a distance of 
about 5" to the SW. A total of 17 excess counts were collected 
(including the area around the faint counter jet) in the energy 
range of 0.4-2.4 keV. 

The X-ray spectrum of DG Tau itself revealed two indepen- 
dent X-ray co mponents in both the Cha ndra and XMM-Newton 
observations dGiidel et all I2005L l2007bb ; a soft, little absorbed 
component is emitted by cool (w3— 4 MK) plasma; and a hard, 
strongly absorbed component originates from hot, occasionally 
flaring (« 20-70 MK) plasma. The soft component was at- 
tributed to emission from the base of the jets. 



2. The target: DG Tau and its jets 

2.1. An accreting, jet-driving T Tauri star 

DG Tau is a single classical T Tau star dLeinert et all Il99lb . 
surro unded by a disk of dus t ( Dutrev et all [19 96: Testi et all 
12002b a nd gas dDutrev et alj|1996t lKitamuraet1d1ll996albl = 
K96a,b; iTesti et al.ll2002b~ The dust disk as detected in the mil- 
limeter continuum is relatively compact (due to limited sensitiv- 
ity of present-day instrumentation to o ptically thin emissi on of 
the outer disk), with a size of 17 1 x 076 dDutrev et al.ll 1996 . also 
K96b). Dust is also present at larger distances as inferred from 
extinction of the counter jet (see below). CO gas remains easier 
to detect at large distances from the star; K96a,b find a 13 CO disk 
size of about 40" x 30" although we note that this structure is 
not a Keplerian disk but more akin to a residual, flattened enve- 
lope. Interestingly, line shifts as a function of position indicate 
that the outer parts of the disk are expanding, perhaps driven by 
interactions between the stellar wind and the disk surface. This 
latter wind would also be responsible for blowing off residual 
envelope gas; a spherical gas envelope is indeed not present any- 
more, as inferred from modeling of the 13 CO line (K96a). 

DG Tau is a strongly accreting CTTS, with an observed ac- 
cretion rate of M - (lO' 7 ^ -J_0_ 613 )M o yr' 1 dWhite & Ghezl 
1200 It IWhite & Hillenbrand. 12004b . It also ejects a well stud- 
ied collimated jet with knots and bow shocks out to at least 



3. New Chandra observations 

We have obtained new observations of DG Tau with Chandra 
ACIS-S (Table[T|i, for a total of ^60 ks of exposure time. We have 
merged these data with the previously obtained Chandra data to 
produce images with an equivalent exposure time of approxi- 
mately 90 ks, i.e., three times as much as reported in Giide l et al.l 
(2005). The new observations were collected in three segments, 
two shorter ones in December 2005 (ObsID 6409 and 7247) and 
a longer one in April 2006 (ObsID 7246). The data reduction 
fo llowed stand a rd Ch andra CIAO analysis thread^] as described 
m iGiidel et alj d2005b . We used the "Very Faint" mode to effi- 
ciently reduce background radiation. DG Tau's optical sky coor- 
dinates for Epoch 2005.0 are RA( 2000.0) = 04h 27m 04 .697s, 
dec(2000.0) = +26deg 06' 16710 dDucourant etalll2005b . 

We extracted stellar X-ray spectra individually for each ob- 
servation using counts within 174 - 177 (depending on the shape 
of the stellar image) of the centroid position of the DG Tau stel- 
lar source. This extraction radius contains approximately 95% 
of the power in the PSF (see Chandra Proposers' Observatory 
Guide [POG] v.9)0 Counts from regions attributed to the jets 



1 http://cxc.harvard.edu/ciao/guides/acis_data.html; CIAO version 
3.3.0.1 was used 

2 http://cxc.harvard.edu/proposer/POG/html/ 
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Fig. 1. Chandra ACIS-S images of DG Tau and its jets. North is up and east is to the left. From left to right: a) Count image for the 
0.6-1.7 keV range; b) same but smoothed; c) color-coded, smoothed hardness image; d) same but with different color coding (see 
text for details). Pixel size is 07492. 



were extracted from circles with radii of 1.75" and 2.6", off- 
set by a few arcs econds to the NE and SW, respectively (see 
Gud el et al.ll2005l for an illustration). 

Spectra were produced for DG Tau, the SW (forward) jet, 
and the NE (counter) jet for each observation (using the CIAO 
specextract task). We fitted the spectra with simple thermal, 
collisional-ionization equilibrium models (vape c) in the XSPEC 
software package (version 1 1.3.1; lArnaudll996l) . complemented 
with photoelectric absorption models. DG Tau required two ther- 
mal models that were each subject to different amounts of ab- 
sorption. The two jet sources were fitted jointly (see below). 
We fixed all element abundances at values that are typical for 
T Tauri stars, as used in the XEST project (Gu del et al.ll2007at 
iTelleschi et al.ll20~07al and references therein)0 

Because the ACIS-S background level is negligible for our 
small source s (see below), we took advantage of the C statistic 
dCashlll979l) for the spectral fits in XSPEC, applied to unbinned 
spectra in the energy range of 0.2-7 keV for DG Tau and 0.2- 
2.2 keV for the jet sources. 



4. The X-ray jets of DG Tau 

Our X-ray images (Fig.Q]) provide clear evidence for both a for- 
ward (to the SW) and a counter-jet (to the NE), symmetrically 
arranged outside the stellar point-spread function (PSF) out to 
distances of about 5". This is the first and so far only young- 
stellar double X-ray jet that can be followed essentially down to 
the star. We now discuss the spatial and spectral results in detail. 



4. 1 . Jet morphology 

An X-ray image of the DG Tau environment is shown in Fig.QJ, 
with a pixel size of 07492. This image was produced by combin- 
ing counts from all four Chandra observati ons. The jet proper 
motio n amounts to about 0715-073 pe r year dEisloffel & M undt 
Il998l for the 1983-1990 time interval: iDougados et al.ll2000l for 
1994-1997). If the proper motion was similar during the time 
interval of our Chandra observations (2004-2006), then the jet 
features would shift by about one ACIS-S pixel between the ear- 
liest and the latest observation; such a shift is not critical given 
the size of the jet source (see below). To minimize the exten- 
sion of the combined image of the DG Tau stellar PSF due to 
slight systematic offsets between the attitude solutions of the 



3 The adopted abun dances are, with respect to the solar photospheric 
abundances given bv I Anders & Grevessd Il989l) : C = 0.45, N = 0.788, 
O = 0.426, Ne = 0.832, Mg = 0.263, Al = 0.5, Si = 0.309, S = 0.417, 
At = 0.55, Ca = 0.195, Fe = 0.195, Ni =0.195. 



observations, we determined the stellar centroid coordinates in 
each exposure using the CIAO wavdetect task, and then shifted 
the centroids to a common coordinate (the maximum shift ap- 
plied to one of the exposures was 0742, or less than one pixel). 
Also, the standard "pixel randomization" procedure was turned 
off (but was left on for the data sets used for spectral extraction). 
Because the boresight coordinates are slightly different for the 
different exposures, the data were reprojected to a common tan- 
gent point, using the reproject_events task in CIAO. 

To suppress background and to emphasize the soft jet 
sources, only counts falling within the 0.6-1.7 keV range are 
plotted in the figure. There is clear evidence for a jet-like ex- 
tension to the SW along a position angle of ^225 deg, but 
we also find a significant excess of counts in the NE direc- 
tion (PA »45 deg). This is coincident with the j et optical axis, 
which for the S W jet has been given as 222 deg dLavallev et all 
119971) . 226 deg dSolf & Bohnil 19931). or. de pending on individ- 
ual knots, 217-237 deg dEisloffel & Mundtl Il998l) . We are not 
aware of any background sources in this region that could pro- 
duce additional X-ray emission. 

Only very few counts from the jet regions fall outside the 
~0.6-1.7 keV energy range; we therefore use this range for our 
statistic s (somewhat differe nt from the energy range 0.4-2.4 keV 
used bv lGiidel et al. 2005). The SW jet contains, within the ex- 
traction circle defined here, 7 and 1 1 counts for the 2004 and 
the combined 2005/06 observations, respectively (including one 
count in the SW jet at 0.58 keV in ObsID 6409). For the NE jet, 
the numbers are 4 and 5 counts, respectively. Based on the 2004 
observation, we would have expected to detect, in the 2005/06 
observations, a total of 14 counts in the SW jet and 8 counts in 
the NE jet, but the differences to the actually observed counts 
are within about lcr of the uncertainties due to counting statis- 
tics, and are therefore not significant. The combined Chandra 
exposures thus collected a total of 18 and 9 counts in the SW 
and the NE jet, respectively (the nearest counts outside this en- 
ergy range were: one count at 0.43 keV in the SW jet of ObsID 
4487; and one count at 2. 1 keV in the S W jet of ObsID 6409). 

We then estimated potential contamination from the diffuse 
(sky and detector) background and from the PSF of the bright 
stellar DG Tau X-ray source. The diffuse background was esti- 
mated by extracting counts from a large area in the vicinity of the 
DG Tau system and scaling them to the jet extraction areas. We 
found the background contribution to be, statistically, 0.6 and 
0.3 cts for the SW and the NE jet area, respectively. We modeled 
the DG Tau PSF using the MARX softwareQ Boresight coordi- 
nates and roll angles were identical to those of the observations, 



http://space.mit.edu/CXC/MARX: version 4.2.0 was used 
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Table 2. Results from spectral fit to jet sources. Errors give lcr 
confidence ranges, and EM and Lx values are for one jet each 



Fig. 2. Comparison of two simulations of the DG Tau X-ray 
sources. Left: The two jet sources are extended; the NEjet source 
is represented by a Gaussian source with sigma = (X'8; the SW 
jet source consists of a Gaussian with sigma =1" plus a linear 
source at a position angle of 225 deg between star and Gaussian 
source. The number of counts in each source is similar to the true 
observations for an exposure time of 90 ks. Only counts in the 
range of 0.6-1.7 keV are shown. - Right: Similar, but both jet 
sources have been defined as point sources. Pixel size is 0'.'492. 



and the simulated source was put at the sky coordinates of DG 
Tau (the simulations were done separately for the 2004 and the 
2005/06 observations, but the results agreed with each other). 
The simulation used the best-fit spectrum of DG Tau except that 
the flux (or the exposure time) was much higher. We then ex- 
tracted counts in the PSF wings using the identical extraction 
regions used for the jets, in the 0.6-1.7 keV range, and found a 
statistical contribution of 0.9 and 0.7 cts of the PSF wings to the 
SW and NEjet source, respectively. We conclude that one count 
per jet area is likely to be due to contamination. This is consid- 
erably less than the Poisson uncertainty in our count numbers. 
These X-ray jet so urces are highly significant (see discussion in 
iGudel et alJl2005l for the first exposure only). A linear feature 
pointing from the stellar PSF to the south, however, appears to 
be due to a coincidental arrangement of only four counts. 

To characterize the jet morphology further, we smoothed the 
image using CIAO task aconvolve, treating the jets and the stel- 
lar PSF separately. In Fig.QJ) we used a Gaussian sigma of 0.75 
pixels for the area outside the star, and 0.5 pixels for the stellar 
image. The smoothed image clearly shows the extended mor- 
phology of the two jet s ources, reminiscent of the optical image 
at least for the SW jet dDougados et al.Ll2000T) . 

Figure [2] shows simulations of two different jet source 
shapes, based on the MARX software. In the left figure, the SW 
and NE jets are composed of a Gaussian with sigma = 1", and 
the SW jet additionally contains a 1-D linear source of 2" length, 
positioned radially outside the stellar PSF at a position angle of 
225 deg. The linear feature contains k 1.4 times as many counts 
as the Gaussian source. The right figure shows a simulation in 
which both jet sources are point sources. 

The total number of counts in each simulated jet feature is 
similar to the numbers in the real 90 ks observation. The ob- 



5 The half-energy radius of the point response function of an 
ACIS-S on-axis source is about d'41 or ss 1 pixel, and the 
80% encircled-power radius is w (X'7 or nearly 2 pixels (see 
Chandra Proposers' Observatory Guide v.9) . We also simulat ed the 
PSF i n the Chandra Ray Trac er (ChART; ICarter et alj [2002) soft- 
ware (http://cxc.harvard.edu/soft/ChaRT/cgi-bin/www-saosac.cgi ) for 
the position of DG Tau and the boresight parameters of ObsID 4487 and 
7246. The PSF was found to be compact, its core being very slightly 
elongated along a SE-NW axis, as can also be seen in Figs. [TJ and [2] 
There is no extension in the directions of the jets (SW-NE). 



SW Jet: 



NE Jet: 



M a J> 

AA'h 
EM* 

L x [0.1-10 keV] 



5.7xl0 21 crrr 



3.4+" x 10 6 K 



-0.8 

1.2 x 10 28 



erg s 



The two Nh values are coupled through the fit parameter AA'h- 
Therefore, no independent errors are given for the NEjet. 
Afiyorward and EM are, within lcr, poorly constrained to low values 
and T is poorly constrained to high values because there is a second 
(probably unphysical) minimum in the statistic, with Miforward ~ 0, 
T = 6.8 MK, and very small EM. See text for details. 



servation resembles the simulation of extended sources; it is not 
compatible with one point source per jet. A more quantitative 
comparison will need longer exposures of these features. 

We next generated "hardness" images, using two comple- 
mentary procedures. In the first case, each individual count was 
attributed a color as a function of its energy. Energies in the 
0.55-0.7 keV range are represented in red, energies in the 1.5- 
1.75 keV range in blue. In the 0.7-1.5 keV range, colors change 
continuously from red to yellow to blue. The resulting image 
was smoothed, using a Gaussian with sigma = 0.85 pixels - see 
Fig. [It. 

For an alternative representation, we extracted three pixel 
images smoothed as before (Gaussian smoothing with sigma = 
0.75 pixels), a "red", "green", and "blue" image for the 0.55- 
0.85 keV, 0.85-1.2 keV, and 1.2-1.75 keV ranges, respectively. 
The combined false-color image is shown in Fig.QJl. 

Both hardness images show that the counter jet is harder, 
with photon energies mostly above 1 keV, while the forward jet 
shows a mixture of softer and harder counts. 



4.2. Spectral analysis of the jet sources 

Only simplistic models can be fitted to the spectra of the jets, 
given the small number of counts in the two jet sources. We ex- 
perimented with isothermal plasma components (using the vapec 
model in XSPEC, with the abundances mentioned above) subject 
to photoelectric absorption (using the wabs model in XSPEC). 
Some plausible assumptions were made to keep the number of fit 
parameters low: Because the jet morphologies and the lengths of 
the two jets are similar, we adopted equal emission measures and 
(for similar heating mechanisms, e.g., shocks with similar veloc- 
ity) equal temperatures. The absorption column densities, Nh, 
remained different for each source because different gas com- 
ponents (e.g., the star's extended circumstellar disk) may absorb 
the jets differently, as also suggested by the hardness images dis- 
cussed above. The initial fit parameters were therefore the emis- 
sion measure (per jet), a single electron temperature (for both 
jets), and two Nh values. We found, however, that the absolute 
values of A^h are poorly constrained due to numerical cross-talk 
with the temperature (lower A^h can be compensated by higher 
temperatures), while the difference ANh = A^counter _ ^.forward 
between the absorption of the the counter jet and the forward jet 
is sufficiently well confined. This is illustrated in Fig. [3] where 
the confidence regions are shown on the A^counter vs. A^forward 
plane. The region of optimum A^h values shows a strong and 
nearly linear correlation between the two values, in such a way 
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Fig. 3. Confidence regions for A^h of the counter jet versus A^h 
of the forward jet, using an isothermal plasma component per 
jet (with equal EMs and temperatures). The contours show the 
lcr (68%, in red), 90%, (in green), and 99% (in blue) confidence 
regions for two interesting parameters based on the C statistic. 
The black contour is plotted to outline the regions of lowest C 
statistic (referring to AC = 1). The cross marks the location of the 
lowest C statistic, although the minimum at A^tbrwai-d = 0.3 and 
A^H.counter = 0.57 is more likely to represent a physical solution. 



that their difference is well constrained. We have therefore cho- 
sen to fit A^H.forward and ANh, while the value of A^counter can be 
derived from any pair of these two parameters. 

The best-fit AA^e is 2.7 x 10 21 cm -2 , with a lcr range of 
(1.6-4.2)xl0 21 cm" 2 and a90% range of (0.9-5.2)xl0 21 cm" 2 . 
A differential absorption column between the two jets is thus 
very likely, and is in fact expected as discussed below. The plot 
shows two minima, the slightly deeper one located at A^h forward ~ 
0.02 x 10 22 cm" 2 and A^counter ~ 0.29 x 10 22 cm" 2 , and a 
slightly shallower one at A^forward ~ 0.30 x 10 22 cm" 2 and 
A^H.counter ~ 0.57 x 10 22 cm" 2 . The two solutions agree at the 
lcr level; however, the former solution requires a high best-fit 
temperature (0.59 keV) which, together with an Ng value that is 
unusua lly low for Taurus pre-main sequence stars (Gude fet all 
2007a), suggests that the solution is unphysical. We adopt the 
second solution as the more reasonable "best-fit" in the follow- 
ing although we repeat that only AA^h can be sufficiently well 
constrained. 

Despite the large errors of the fit parameters, three features 
are noteworthy. First, the absorption column density toward the 
forward (SW) jet is small (« 3 x 10 21 cm" 2 ), agreeing with A^h of 
the soft but not the hard stellar component (see below). Second, 
the absorption toward the counter jet is higher, compatible with 
the increased hardness discussed in Sect. 14.11 And third, the elec- 
tron temperature of the jet sources is low («3.4 MK) compared 
with c oronal temperatures of T Tauri stars (e.g.. iTelleschi et all 
I2007bl) . 

We do not give errors for Lx as these depend very sensitively 
but non-tr ivially on the rather uncertain absorption column den- 
sities (see Giidel et al. 2007a for estimates of lower limits to the 
uncertainties of Lx of many X-ray sources in the Taurus star- 
forming region). Based on the best-fit values, the total X-ray 
output from the jets outside the stellar PSF (radius of 1'.'4 - 177) 
would amount to a few times 10 28 erg s" 1 or 10-20% of the stel- 
lar soft component (see below). 



J fa | 

t+ V/u dgtgu 




0.5 



1 2 
channel energy (keV) 



Fig. 4. Co-added ACIS-S spectra of the DG Tau point source 
(upper spectrum, black; the lowest-energy bins are connected 
for clarity) and the two jets (lower spectra; spectrum of less ab- 
sorbed forward jet is shown solid, in red, while spectrum of the 
more strongly absorbed counter jet is shown dotted, in blue). For 
the jet sources, the solid histogram shows the joint spectral fit. 
The spectrum of DG Tau has been binned such that at least 15 
counts are contained in each bin; the jet spectra have been binned 
arbitrarily to produce non-zero bins from ss 0.3 keV to 2 keV. 

5. The DG Tau X-ray source 

The two-component spectral phenomenology of TAX sources 
dGiidel et all 12007b) is present in each of the new Chandra ob- 
servations. The spectral-fit results are reported in Table [3] to- 
geth er with the re s ults fro m the XMM-Newton observation taken 
from iGiidel et all d2007b). The rather short exposures obtained 
in 2005-2006 led to considerable uncertainties in the derived pa- 
rameters. This is particularly evident in the absorption column 
density, N-a,i, the electron temperature, kT\, and consequently 
the X-ray luminosity, Lx,i> of the softer spectral component for 
the shortest observation, ObsID 7247, comprising only 142 cts. 
On the other hand, the count rates of the soft component agree 
within the error bars for all Chandra observations (Table O, and 
no variability was recorded in the soft component during the in- 
dividual exposures. Because the 90% error bars of Nn,i and kT\ 
also strongly overlap for the four Chandra exposures, we per- 
formed a joint fit of these spectra in the energy range of 0.2- 
1.1 keV, reported in the penultimate column in the table. The 
resulting parameters of the soft component compare very favor- 
ably with the XMM-Newton EPIC PN results (last column). 

Because the hard component is variable (Table [3j slow, non- 
periodic modulations by a factor of ^2 were present in ObsID 
7246 although larger flares were absent), agreement between 
the fit parameters is not expected, but we notice that the hydro- 
gen absorption column density, A^h,2, agrees in all observations 
within the 90% error ranges. It is remarkable that A^h,2 is approx- 
imately 20 times higher than A^h.i- 

For illustration purposes, we show the combined Chandra 
spectrum in Fig. [4] but again emphasize that spectral fits were 
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Table 3. Results from spectral fits for the DG Tau point source. Errors give 90% confidence limits 



Chandra: XMM-Newton: 



KJUSLLf 


AA9.1 




77/L7 


77A^ 


3.11 Chandra 


rnrn^zLrnni 

UZUJJ^uZUl 




Tan 9004 


rw 7005 




Anr 7006 


9004-9006 
zuuu 


Aim 7004. 


counts [0.2-7 keV] 


388 


204 


142 


349 




634 


count rate [ct s -1 , 0.2-7 keV]* 


0131(7) 


o oi26(9i 


0089(7) 


0.0125(7) 




0206(10) 


Soft component: 














count rate [ct s~\ 0.2-1.1 keV]* 


0.0046(4) 


0.0041(5) 


0.0041(5) 


0.0046(4) 




0.0083(6) 


AT Hjl [10 21 cm- 2 ] 


1 4+1.8 


95 +2 3 


4 5 +41 


2.0 +2 « 


1 3 +u 

1 -'-0.6 


\.\* lA 


T } [MK] 


4 6 +2 3 

^'-13 


6.0 +1 7 


3 2 +15 


3 6 +2:3 


4 5 +0 - 7 

^-0.9 


Q 7+81 
J - '-0.5 


L x ,i [ergs" 1 , 0.1-10 keV] 


1.1 x 10 29 


8.2 x 10 28 


4.5 x 10 29 


1.7 x 10 29 


1.0 x 10 29 


0.96 x 10 29 


Hard component: 














count rate [ct s~\ 1.7-7 keV]' 


0.0063(5) 


0.0068(6) 


0.0031(4) 


0.0066(5) 




0.0099(7) 


7V H ,2 [10 21 cm" 2 ] 


21+ 6 


^9+16 

->^-n 


15 +14 


33 +n 

-u 




181 


T 2 [MK] 


31+ 21 


27+25 


34 +77 


23+23 




6 9+ii9 

u7 -25 


L x ,2 [ergs- 1 , 0.1-10 keV] 


6.9 x 10 29 


1.0 x 10 30 


2.8 x 10 29 


1.1 x 10 30 




9.6 x 10 29 



a Average properties for the hard component during XMM-Newton EPIC PN observations (including flares); counts for 24932 s of low- 
background data 
b parentheses give errors in last digit 



performed for the individual, unbinned s pectra (illustrations of 
spectral fits to binned data are sho wn in iGiidel et al.l 120051 for 
Chandra and in lGiidel et aT]|2007bl for XMM-Newton). 



6. Discussion 

Our Chandra observations add unprecedented information to the 
X-ray emission model of DG Tau. In particular, we report the 
detection of an extended, bipolar X-ray jet down to the stellar 
PSF and out to a distance of about 5". This is the first double- 
sided X-ray jet reported from a pre-main sequence star. 

The ACIS-S images reveal the presence of jets in the SW and 
NE directions, coincident with the direction of the optical jets. 
The NE jets appears to be slightly harder and also shows fewer 
photons in a similar extraction region. Further, the spectra of the 
DG Tau point source reveal a TAX spectral morphology. The 
DG Tau system thus hosts at least four distinct X-ray sources of 
different origin and subject to different gas absorption columns, 
as illustrated in the sketch shown in Fig. [5] These components 
are: 

1. a weakly absorbed, diffuse, soft component along the 
forward-jet axis; 

2. an more strongly absorbed, diffuse, soft component along 
the counter-jet axis; 

3. a weakly absorbed, compact, non-variable, soft component 
(within the stellar PSF); 

4. a strongly absorbed, compact, flaring, hard component 
(within the stellar PSF). 



We now discuss the various features seen in the X-ray spectra 
and images, and propose a model consistent with all observed 
features. 



6.1. The X-ray jets: Physical properties 

In this section, we will be interested as to how the jet cools, 
and why no X-ray emission is seen at distances beyond 5" from 
the star. There are further prominent bow shocks outside the re- 
gions discussed here, in particular a structure that was located at 
l(y/7 from DG Tau in ob servations obtained in 1986 December 
dEisloffel & Mundtt 1 19981) and still within a « 12" a decade 
later dStapelfeldtl 1 19971) Using a proper motion of 07 15 -07 3 
per year (Eisl offel & Mundtl ll998HDougados et all |2000|) . this 
bow shock should now have expanded to about 15". Obviously, 
the DG Tau jet consi sts of a series of knots th at are frequently 
ejected from the star {Eisloffel & Mundtl [T998I) . If the more dis- 
tant knots were similar to those now seen in X-rays when at simi- 
lar distances from the star, then the heated gas must have cooled, 
or the emission measure have decreased below our detection 
limit because we have found no significant X-ray source at those 
distances. Given the time to travel from » 5" to « ll"(the re- 
gion behind the outer bow shock, in the "X-ray free" lower right 
corner of Fig. [TJ, we adopt an X-ray decay time (with respect to 
the range of sensitivity of the X-ray detectors) of « 20 - 40 yrs. 

To explore the potential cooling mechanism(s), we first dis- 
cuss two extreme cases of cooling: cooling by expansion, i.e., 
adiabatic cooling; and radiative cooling without expansion. We 
will then discuss cooling if contributions from both cooling 
mechanisms contribute. 

6.1.1. Cooling by expansion 

The DG Tau jet significantly expands with increasing distance 
from the star. Transverse expansion has been measured; vari- 
ous features within a few arcsec of the star show an unusually 
large opening angle in the range of /3 ^ 11-27 degrees (b ased 
on fu ll width at half maximu m, after iDougados et al]|2000l see 
also [Eisloffel & Mundtlll998l) . In contrast, radial expansion (i.e., 
stretching of volumes in radial direction, along the linear flow) is 
unlikely to be important as this would require strong jet acceler- 
ation at distances of several arcseconds, or selective deceleration 
of regions closer to the star, for which no evidence has been re- 
ported. 
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Fig. 5. Sketch of the proposed model of the DG Tau environ- 
ment. Features are not drawn to scale. The observer is located 
toward the upper left. Four X-ray source regions are shown in 
red (the stellar coronal source, and asterisks symbolically mark- 
ing emission regions in the jet) and their absorbing media are 
schematically shown. X-rays seen by the observer are sketched 
by the wavy lines, described on the left, from top to bottom: X- 
rays from the weakly absorbed, spatially resolved X-ray source 
of the forward jet; from the spatially unresolved but spectrally 
identified soft X-ray source closer to the base of the jet; from 
the hard, stellar coronal emission absorbed by the infalling ac- 
cretion streams; and from the spatially resolved X-ray emission 
from the counter jet, slightly absorbed by the intervening disk 
gas. The dashed circle denotes the observing PSF of the star; 
features inside the PSF are unresolved from the stellar X-ray 
source. 



There are two consequences of transverse expansion: a de- 
crease of the EM and adiabatic cooling (in the absence of radia- 
tion), both leading to a decay of the X-ray emission in time. 

The cross-section area of the jet at distance r from the star is 



or for small changes in A(r), 
dA dr 

T = 2 T- 



(2) 



For the time span of 2 years, dr = (X'3-0'.'6, and therefore 
dA/A kQ. 12-0.24 at a distance of 5". The electron density de- 
creases by the same factor as the volume increases, and therefore 
the emission measure, EM = n^V, decreases by 12-24% at a dis- 
tance of «5" within 2 years. Such a decrease is too small to be 
significantly measured in our 2004-2006 data. 

However, for a distance of, say, 8"-l 1" (the lower-right area 
in FigQ}, A(r)/A(ro) - (r/ro) 2 ~ 2.6^1.8. The emission measure 
for a given mass element in the jet thus decreases by a factor 
of 2.6-4.8 from a distance of 5" to 8-11" from the star due to 
expansion alone. Because the gas is optically thin to X-ray emis- 
sion, the X-ray surface brightness decreases by a factor of 2.6 2 - 
4.8 2 = 6.8-23, making its detection at this distance against the 
background level very difficult in our exposures. 

At the same time, the plasma cools by expansion. In the lim- 
iting case, we neglect radiative losses, i.e., the gas cools adiabat- 
ically. Then, 

TV 7 ' 1 = constant (3) 
with y — 5/3 for monatomic gas, and for small changes in T, 
2dA Adr 



T V 3 A 3 r 



(4) 



Using dA/A « 0.12 - 0.24 for a two-year interval, we find 
\dT/T\ ^ 0.08 - 0.16. Again, this temperature decrease is too 
small to be measured during our observing interval. 
For long intervals, however, 



m = / aw r 2/3 = (rv 4/3 

T(r ) \A(r )) \r ) 



(5) 



7 9 P 

A(r) = nr tan -. 

2 



(1) 



which, for r =8"-l 1", the above projected jet velocities and ro = 
5" yields T(r)/T(r ) ~ 0.35-0.53, i.e. T(r) *1.2-1.8 MK. 

A temperature decrease will also affect the count statistics 
in the detector, because the effective area decreases toward the 
softer energy range. We simulated the specttum of the forward 
jet in XSPEC, adopting a much higher flux for better statis- 
tics and further changing its temperature. We then measured the 
simulated count rate in the 0.6-1.7 keV range. Lowering kT 
to 0.249-0.272 keV (T = 2.9-3.2 MK), which corresponds to 
\dT/T\ =s 0.08-0. 16, reduces the count rate by 1 1-26%. Together 
with the count rate drop from the decrease of the EM (see above), 
the total drop of the detected count rate would thus be » 23- 
50%. It is possible that this decrease affected our measured count 
rates, but we cannot prove such cooling effects as the count-rate 
decrease was not statistically significant during the two years 
covered by our observations. Furthermore, as some of the ma- 
terial in the detected jet structure may cool, new hot material is 
likely to be replenished from closer to the star, as we have found 
the X-ray structure to be extended toward the stellar PSF. 

On the other hand, considering a longer time span of 20 
years, i.e., and expansion from 5" to 8"-ll", therefore lower- 
ing the plasma temperature to 1.8 MK, the count rate decreases 
by 83%. This, together with the sttong decrease of the EM, re- 
inforces our finding that a cooling, expanding X-ray jet starting 
under presently observed conditions at 5" distance will be diffi- 
cult to detect at distance s of 8-1 1", i.e., the regio n just behind a 
bow shock described bv lEisloffel & Mundtl d!998l) . 

In summary, adiabatic expansion alone should lead to the 
disappearance of the X-ray jet at distances significantly beyond 
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5" in Chandra exposures such as ours. Our exposures are not 
sufficient, on the other hand, to verify a cooling effect during the 
two-year time span. 

6.1.2. Cooling by radiation 

Radiative cooling may be significant as well. We cannot esti- 
mate the radiative cooling time because the plasma density is 
unknown, or in other words, it is not clear whether the observed 
X-ray radiation extracts a significant fraction of the thermal en- 
ergy of the hot plasma within the 2-20 yrs of interest here. We 
will, on the other hand, study consequences for the X-ray jet 
radiative cooling is the dominant loss mechanism. The volume 
of the X-ray detected source of the forward jet is approximately 
(see Sect, gj} 



(6) 



3" * (075) 2 tt / „\3 , 

sin(38 deg) v ' 

where 38 deg is the jet inclination (Eislo ffel & MundtL fl998). 
we have adopted a length outside the PSF and a cross-sectional 
radius of the X-ray jet of 3" and 075, respectively (Sect. I4.ll ). 
and the final constant gives the conversion from arcseconds to 
cm at the distance of Taurus. We thus find V = 3.6 x 10 46 cm -3 
but note that the X-ray gas may occupy only a fraction of this 
volume, defined by the volume filling factor /. 

With the best-fit EM from Table [2] the density of the X-ray 
emitting jet gas is 



'EM 

Jv 



1/2 



170 



1/2 



[cm" 3 ]. 



(7) 



The thermal energy decay time then is, assuming that energy is 
lost by radiative cooling (i.e., decrease of T) only, 



T = 



3kT 
n~MT) 



* 1.4 x 10 3 r/ 1/2 [s] * 15400/ 1/z [yr] 



1/2 



(8) 



where A(T) a 1.7 x 10~ 23 erg cm 3 s _1 is the cooling func- 
tion (evaluated at 3.4 MK for our model, derived from apec for 
the 0.001-100 keV interval), and in the last equation we have 
adopted the X-ray measured jet temperature of 3.4 MK. 

Here, the energy decay reflects exclusively in a decrease of 
the temperature, T, while the emission measure remains con- 
stant. As before, we model the count rate reduction for decreas- 
ing temperatures for the ACIS-S detector. Given that we detected 
18 cts in the forward jet, we now seek the temperature for which 
the same emission measure results in a detection limit of only 5— 
6 detected counts, i.e., a count rate reduced to 25-30%. For fewer 
counts, the source will be difficult to detect, depending also on 
the spatial distribution of the counts (we recall that approxi- 
mately one count will be due to contamination from background 
radiation and the stellar PSF). Using the detector response, we 
found that a count rate reduction to 25-30% corresponds to a 
temperature reduction to 50-61% of the initial value, i.e., a time 
interval of 0.49-0.56 e-folding decay times of the energy (Eq.[8]). 

If the jet should become undetectable in an observation like 
ours after 20-40 yr, i.e., for an energy e-folding decay time of 
10-22 yr, we require / = 4 x 10~ 7 -2 x 10~ 6 , or densities of 
(1.2 - 2.6) x 10 s cm" 3 . 

6.1.3. Cooling by expansion and radiation 

If radiative losses are significant during the time span of inter- 
est here, then the adiabatic approximation in Sect. I6.1.T1 breaks 
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Fig. 6. Cooling behavior of the jet plasma, assuming transverse 
expansion as observed in the optical jet, and radiative losses for 
various initial electron densities, no = n e (t = 0) at a stellar dis- 
tance of ro = 5". The curves are color-coded for initial den- 
sities given in the plot. The two curves for no — 10 cm -3 and 
no = 100 cm -3 are nearly coincident. 

down. On the other hand, if we accept that the hot plasma is sub- 
ject to expansion as observed for the cooler optical jet, then the 
treatment of radiative losses in Sect. 16.1.21 is also not sufficient. 
Both cooling terms must be combined. 

We consider the first law of thermodynamics, 



dU + SW = 6Q. 



(9) 



Here, dU is the change of the internal energy, U = 2aN e kT, 
where N e is the total number of electrons and a = 1 / (y— 1 ) = 3/2 
for monatomic gas (for simplicity, we adopt a hydrogen plasma 
with a total number of particles N - N p + N s - 2N e ; the differ- 
ence to a realistic plasma is not significant, the electron density 
in the latter being only about 10% larger than the ion density). 
Further, 6W = pdV is the pressure work done by the gas, with a 
pressure of p — 2n e kT . If radiative losses are present, we set 



SQ = -n e (t) 2 V(t)A(T)dt 



(10) 



where A(T) is the temperature-dependent cooling function (in 
erg cm 3 s" 1 ) and n e is the electron density. Substituting these 
expressions into Eq. (O, we find 



dT dV 

a W) + W) 



n e (t)A(T) 



dt. 



(ID 



2kT(f) 

Note that in the absence of radiative losses, Eq. (fTTT i is equivalent 
to equations [3] and |4] for an adiabatic process. Radiation makes 
the problem explicitly time-dependent; the initial electron den- 
sity, no = n e (f = 0), is introduced as a free parameter. 

We define the observed expansion of a given volume element 
by the radial distance along the jet axis, r(f) = ro + vt. Then, 
V 



r 2 

oc y- { 



cv /- (in the absence of expansion along the jet axis), and 



»-2 



Substituting these expressions, we find 



a- 



dT 

W) 



2v IroY n Q A(T) 
r(t) + \r(t) 2kT{t) 



dt. 



(12) 



We have determined the run of T numerically, starting with 
To = 3.4 MK and with various no values, taken to be the ini- 
tial electron density at ro =5" (not yet considering a strict lower 
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limit to «o - see below). The cooling function was approximated 
by A(T) = 4.83 x 1(T 23 - 9.75 x 1(T 24 T erg cm 3 s 4 (T in 
MK) for the range of T =1-3.5 MK, as found from an in- 
tegration of the spectral energy distribution for an isothermal 
plasma in XSPEC, using the vapec code with our adopted el- 
ement abundances (Sect. [3]). For T < 1 MK, we kept A constant 
at 4 x 10~ 23 erg cm 3 s . We used an average projected jet ve- 
locity of 0.225" yr _1 . The results are shown in Fig. [6] We see 
that - in agreement with the estimates for pure radiative losses - 
the cooling process is dominated by expansion up to densities of 
several times 10 3 cirT 3 , i.e., cooling is adiabatic. Only for den- 
sities S 10 4 crrT 3 does radiation matter on time scales of tens of 
years. 



6.1 .4. Implications from jet cooling 

Because the electron density in the jet X-ray source is unknown, 
the importance of radiative cooling cannot be directly assessed. 
We will, however, discuss the limiting case of low densities. 
Equation (Q gives the lowest density for the maximum volume 
filling factor, / = 1, namely n e = 170 cm 4 . This implies an 
upper limit for the thermal energy of the hot gas of 3n e kTV ~ 
8.6 x 10 39 erg in the forward jet alone or « 1.7 x 10 40 erg in both 
jets together. 

The mass outflow rate i n the DG Tau jets is 6 X 10~ 8 M yr _1 
(twice the value given by Lavallev-Fouqu et et al.ll2000l for the 
forw ard jet) with a characteristic flow velocity of 30 km s _1 
(e.g., lEisloffel&Mundllll998t iDougados et all 120001) . The ki- 
netic energy rate is therefore Aifj e tV? t /2 « 1.7 x 10 33 erg s _1 . 
The total jet X-ray luminosity is, according to Table [2j 2.4 X 
10 28 erg s _1 , i.e., a fraction of 1.4 x 10~ 5 of the kinetic energy 
dissipates in X-rays after transformation to thermal energy of the 
hot plasma. 

The jet out to a distance of 5" forms within m 20 yrs based 
on the measured proper motion. The mass outflow rate then in- 
dicates that ss 1.2x 10~ 6 M are contained in the jets, with a total 
kinetic energy of 1.1 x 10 42 erg. Therefore, a maximum of 1.6% 
of the kinetic energy is transformed into thermal energy in the 
hot plasma. 

We can now investigate whether the pressure in the hot gas 
contributes to the expansion of the jet. The bulk gas observed in 
optical lines has a density of 10 3 cirT 3 - 10 4 cm -3 at distances of 
a few arcs ec from the star, with an ioniz ation fraction approach- 
ing unity (iLavallev-Fouauet et all l2000b . The temperature is of 
order 10 4 K (Hamann, 11994 . The gas pressure is therefore of 
order P bu i k * (3 x 10~ 9 - 3 x 10~ 8 ) dyn cm 4 . For the X-ray 
emitting gas, we estimate, from the above values, a minimum 
pressure for / = 1, namely Phot ~ (1 -6 x 10~ 7 ) dyn cirT 2 . A fill- 
ing factor of unity is impossible because only a small fraction of 
the energy heats the hot plasma while cooler gas predominates. 
It thus appears that the hot gas pressure is grossly out of equi- 
librium with its environment, contributing to the transverse jet 
expansion discussed in this section. 

We emphasize that this argumentation holds regardless of the 
exact filling factor or of the relevance of radiative cooling. It is a 
minimum-density estimate related to the energy decay in Eq. ([8j. 



6.2. Heating the X-ray jets 

How are the jets heated to the observed temperatures? Most of 
the conventionally detected emission from jets of pre-main se- 
quence stars comes from low-ionization transitions such as [O i], 
[Nn], or [Sn], indicating temperatures of no more than a few 



thousand K. Hotter gas has been identified, e.g., in emission 
from He i, \Q ml , C m, and O iv, giving evidence for hot winds 
(ICohen & Fullerlll985tlTakami et all[2002l;lDupree et al.Ll2005l) 
although some of these claims have been questioned, and max- 
imum wind temperature s of on l y 2 x 10 4 K have b een derived 
dJohns- Krull & Herczegl 120071; Kwan et all 120071) . Measured 
shock velocities in the je t s of D G Tau are of order 50-100 km s _1 
(lLavallev-Fouquet et all 120001) . sufficient to heat gas to a few 
10 s K but not to the observed 3-5 MK. 

A possibility is that a fraction of the jet gas collides di- 
rectly with the interstellar medium, in which case the shock ve- 
locity would be several 100 km s and therefore in principle 
sufficient to heat a fraction of the gas to X-ray emitting tem- 
peratures. Alternatively, magnetic heating may be in operation. 
Magnetic fields are commonly invoked to explain jet accelera- 
tion near the star and the disk, and also to coUimate jets during 
their propagation dBlandford & PavneLll982tlUchida & Shibatal 
fT985llKonigl & Pudritzi l2000tlShu et all 120001) . 

Ambipolar diffusion heating results from the separation of 
charged particles from neutrals acro ss magnetic fi e lds. The pro- 
cess has been described in detail bv lSafierl d 19921) . iGarcia etail 
(12001 ah . and lGarcia et al. (2001b). It is most efficient as long as 
the gas is weakly ionized and therefore not relevant for our hot 
plasma. Self-consistent jet models converge to jet temperatures 
of « 10 4 K due to ambipolar diffusion heating. 

Ohmic dissipation of currents operates in highly ionized gas. 
Initially jet-aligned magnetic fields are wound up due to rota- 
tion, producing helical fields which drive currents. Depending 
on V x B, sufficient dissipation may be achieved to heat the gas 
to high temperatures. In order for Ohmic dissipation to become 
effective, the gas should be pre-heated, which could be achieved 
by shocks. Direct magnetic-field measurements in the jets of 
DG Tau would be important to quantitatively assess the heat- 
ing mechanism, e.g., based on rad i o synchrotron emission from 
accelerated electrons dCuriel et allll993MRav et"allll997l) . 



6.3. The gas-to-dust ratio in the outer disk of DG Tau 

The hardness difference between the forward and the counter jet 
is easily explained by the presence of a gas disk that absorbs 
the softest photons from the counter jet (see Fig. [5]). This ob- 
servation therefore in principle offers the opportunity to mea- 
sure the gas-to-dust ratio in a circumstellar disk (at distances of 
a few hundred AU) from differential absorption and extinction 
measurements. The spectral-fit results summarized in Table [2] 
show a difference between the gas column densities of the two 
jets of ANh ~ 2.7 x 10 21 cm -2 , albeit with large errors. For 
standard assumptions pertaining to the interstellar medium (e.g., 
dust grain size distribution, atomic gas absorption with "solar" 
eleme ntal composition), Nu =s 2 x 10 21 Ay [cm -2 ] ( Vuong et al., 
2003). Therefore, ANh corresponds to a visual extinction differ- 
ence of AAy 1 .4 mag (with a lo - range of « 0.8-2. 1 mag). The 
difference between t he visual extinctions of the two jets is also 
not precisely known: ILavallev et alj d!997h estimate a difference 
of about 3 mag at a distance of 1" from the star, which would 
suggest nea rly standa rd gas-to-dust ratios when combined with 
AN H , while lPvo etaT] (l2003) propose a difference of as much as 
14.2 mag at the same position. Better determinations of differen- 
tial gas absorption and visual extinction are needed; we also note 
that most of the jet-related X-ray counts were collected from 
larger stellar distances (2"-5") than the distances to which the 
AAy values refer. 
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6.4. The hard stellar source 

Because the hard component occasionally flares, in one case be- 
ing preceded by U b and emission as in solar and stellar flares 
(Giid el et all l2007bl) . it is most straightforwardly interpreted as 
coronal or "magnetospheric". The high temperatures are also not 
consistent with shock heating given gas flow velocities of no 
more than a few 100 km s" 1 (as observed in the jets, or inferred 
from free-fall onto the stellar surface). 

The hard component of the DG Tau point source is unusually 
strongly absorbed, with A^h ~ 2 X 10 22 cm' 2 . F rom the stellar 
extinction, Ay = 1.5-3 mag (Gii del et al.ll2007bl and references 
therein), we expect A^h ~ (3 — 6) X 10 21 cm" 2 . We therefore find 
from Table [3] that the dust content in the accretion streams is 
depleted by factors of 3-6. 

The excess photoelectric absorption requires the presence of 
gas that is depleted of dust in order to suppress optical extinc- 
tion. One possibility are relatively cool stellar or disk winds. 
However, as they expand to large stellar distances, they would 
also affect the jet components and the soft "stellar" component 
unless these were placed very far from the star. Also, it is not 
clear why such winds would be dust-depleted, or if they initially 
are, why they would not form dust. 

Dust is destroyed at several stel lar radii due to stell ar ir- 
radiation; for DG Tau in particular, Giide l et ail (l2007bl) esti- 
mated the dust sublimation radius at (7—10)7?*, similar to the 
corotation radius (see below). The region inside the corota- 
tion radius is thought to be dominated by magnetic accretion. 
The excess absorption is thus most easily explained by the in- 
falling, dust-depleted massive accretion streams. These observa- 
tions therefore provide indirect evidence for dust-depleted accre- 
tion streams that absorb the X-ray emission from the underlying 
corona. 

We now show that the amount of photoelectric absorption 
is plausible. Estimating the gas column density along the line 
of sight through the accretion streams to the stellar corona re- 
quires knowledge of the accretion geometry, which is not avail- 
able. As a conservative limit, we assume spherically symmetric 
infall. The corotation radius, at which the mass begins to accel- 
erate toward the star, can be derived from Kepler's law: 



I GM,P 
. An 2 



2\l/3 



= 9.2xl0 11 cm = 5.6fl», 



(13) 



where M» = O.91M is the stell ar mass dBriceno et all 12002), 
R = 2.46/? is the stell ar radius dGiidel et alll2007ah . P = 6.3 d 
is the rotation period dBouvier et all 1 1 993b . and G = 6.673 x 
10" 8 dyn cm 2 g~ 2 is the constant of gravitation. For a semi- 
circular funnel stream guided by magnetic fields, the highest el- 
evation is reached midway between the corotation radius and 
the stellar surface, at r — 3.3/?„. There, the stream's velocity, 
v; n = 1.3 x 10 7 cm s" 1 if free-fall acceleration is assumed (as 
an upper limit) starting from the corotation radius, is about 40% 
of the fall velocity at the stellar surface, and the ionization de- 
gree of the gas may still be moderate. We note that the jet-axis 
inclination and by infe rence the most likely stell ar rotation-axis 
inclination is «38 deg (Eisloffel & Mundtl TT998). For an order- 
of-magnitude estimate, we assume spherical infall at this radius, 
over a radial distance of 1R*. The mass accretion rate in this case 
is, for spherical symmetry 



M as 47ir 2 «H'« P Vin ~ 8.9 x 10 7 n H 



(14) 



where «h is the hydrogen number density, and m p is the proton 
mass; we assume a mean mass per particle of m p . The observed 



accre t ion rate is M = (10~ 7 - 3 4 - 10~ 613 )M B yr 1 dWhite & Ghezl 
2001; IWhite & Hillenbrand 12004b . We thus find n H ~ 3.3 x 
10 10 cm _3 -5.3xl0 n cm" 3 . Integrated over one stellar radius, the 
absorption column density is Nh ~ 5.7 x 10 21 - 9.1 x 10 22 cm" 2 , 
in agreement with the measured values of order 2 x 10 22 cm" 2 . 

More confined or slower accretion streams will produce 
higher densities, while ionized gas will reduce photoelectric ab- 
sorption. We emphasize that the high accretion rate is crucial to 
obtain high A^h in the above estimates. All TAX sources reported 
before (Gi idel et alll2007bl) are indeed very strongly accreting T 
Tauri stars. 

6.5. The soft "stellar" source 

The soft emission commonly seen in X-ray spectra from other 
T Tauri stars, formed at a few MK as part of a wide distri- 
butio n of plasma in the magnetic corona (e.g., Preibisch et al. 
2005), is not detected here but is absorbed, leaving only a "hard" 
coronal component above 1.5 keV in our spectra. In contrast, 
we see a very strong, separate soft component that is not re- 
lated to the coronal spectrum as judged from its very different 
absorption and the absence of variability during flares seen in 
the co ronal component at higher photon energies dGiidel et all 
2007b). The strong absorption of the coronal X-ray component 
makes an origin of the soft spectral emission from a location 
close to the stellar surface unlikely. In fact, our new observa- 
tions suggest that the gas absorption of the soft X-ray com- 
ponent, Nh = 1-3 (0.7 - 2.4) x 10 21 cm" 2 (90% error range; 
Table [3]l, is lower than the absorption suggested from the visual 
extinction to the star Af H (Ay) «T(3 - 6) x 10 21 cm" 2 (based on 
Ay « 1.5-3 mag, see I White & Ghed200lllWhite & Hillenbrand! 
2004, iMuzerolle et al.lll998l and lHartigan et alj|1995l) . A likely 
origin of these X-rays is the base of the forward jet. Such an 
origin is suggested b y i) the unusually sof t emission not usually 
seen in T Tauri stars dGiidel et all 12007a) . ii) the low Nn,i, and 
iii) the explicit evidence for jets in the Chandra image. Further, 
Nn,i and kT\ agree with the corresponding values of the forward 
jet (Table |3). We therefore suggest that the jets continue to be 
X-ray sources to smaller stellar distances, producing an order 
of magnitude more soft X-rays within the Chandra point-spread 
function of « 2" radius than outside this radius. 

7. Summary and conclusions 

We have unambiguously detected a bipolar X-ray jet associated 
with the strongly accreting classical T Tauri star DG Tau. This 
is the first bipolar X-ray jet reported from a pre-main sequence 
star. The jet is extended at Chandra 's spatial resolution and can 
be followed down into the PSF of DG Tau itself. The jets are 
roughly symmetric as far as length and structure are concerned, 
reaching out to about 5" from the star. 

The X-ray emitting gas of the jet is relatively cool, T « 
3.4 MK, which however still poses a problem for our under- 
standing of the heating mechanism. Shock velocities in this jet 
are too small to heat gas to such temperatures. Ohmic heating by 
magnetic-field driven, dissipating currents may be an alternative. 

We find various gas absorption columns toward the four X- 
ray components detected in the DG Tau system. The soft stel- 
lar component shows A^h smaller than the value derived from 
the stellar visual extinction assuming standard gas-to-dust ratios. 
This suggests that the soft spectral emission originates from a re- 
gion "in front" of the star, i.e., we identify the soft component 
with X-ray emission similar to that from the forward jet but pro- 
duced too close to the star to be resolved in the Chandra images. 
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In contrast, the counter jet suggests stronger absorption, 
which is expected because its X-rays traverse the extended outer 
gas disk. A determination of the gas-to-dust ratio is therefore in 
principle possible by measuring the differential absorption and 
extinction of the two jets. We have succeeded in determining the 
difference in A^h although not the absolute values of each col- 
umn density. The gas column of the intervening extended disk 
structure is relatively small (2.7 x 10 21 cm -2 ), corresponding to 
1.5-2 magnitudes of visual extinction for standard gas-to-dust 
ratios. An extin ction difference of 3 m ag has been reported in a 
previous paper dLavallev etallll997l) . i.e., the outer-disk com- 
position is nearly compatible with interstellar gas-to-dust ratios. 
However, more accurate determination both of the absorption 
difference and the visual extinction difference is needed. 

Finally, the hard component in the stellar spectrum is at- 
tributed to a (flaring) corona. The excess photoelectric absorp- 
tion is ascribed to accretion gas streams that are dust depleted 
because the dust destruction radius is similar to the radius of the 
inner disk edge, thought to be close to the corotation radius. 

Why is the detection of X-ray jets important? The combined 
power of the resolved jets and the unresolved soft spectral com- 
ponent is of order 10 29 erg s _1 or similar to the X-ray output of 
a relatively X-ray faint T Tauri star. This emission is distributed 
above the inner accretion disk. It is therefore an important con- 
tributor to X-ray heatin g and ionization of gaseous disk surfaces 
dGlassgold et all 12004b . a role that has been studied earlier in 
the case of active galac tic nuclei (in the context of the "lamp- 
post model", see, e.g.. lNavakshin &"Ka llman 2001). We specu- 
late that protostellar jets in general develop the same kind of jet 
X-ray emission, but these sources remain undetected close to the 
star because of strong photoelectric absorption. In those cases, 
jet X-rays may act as a dispersed ionization source to affect a 
larger volume of gas than the stellar coronal source alone. 

Acknowledgements. M. A. acknowledges support from a Swiss National 
Science Foundation Professorship (PP002-1 10504). M. A. and S. S. acknowl- 
edge support by NASA through CXC award SAO GO6-7003. The CXC X-ray 
Observatory Center is operated by the Smithsonian Astrophysical Observatory 
for and on behalf of the NASA under contract NAS8-03060. 



References 

Anders, E., & Grevesse, N. 1989. Geochim. Cosmochim. Acta, 53, 197 
Anglada, G. 1995, RevMexAA, 1, 67 

Arnaud, K. A. 1996, in ASP Conf. Ser. 101, Astronomical Data Analysis 
Software and Systems V, eds. G. Jacoby & J. Barnes (San Francisco: ASP), 
17 

Bacciotti, F, Mundt, R., Ray, T. P., et al. 2000, ApJ, 537, L49 
Bacciotti, F., Ray, T. P., Mundt, R., Eisloffel, J., & Solf, J. 2002, ApJ, 576, 222 
Bachiller, R. 1996, ARA&A, 34, 111 
Bally, J., Feigelson, E., & Reipurth, B. 2003, ApJ, 584, 843 
Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 
Bouvier, J., Cabrit, S., Fernandez, M., Martin, E. L., & Matthews, J. M. 1993, 
A&A, 272, 176 

Briceno, C, Luhman, K. L., Hartmann, L., StaufFer, J. R., Kirkpatrick, J. D. 

2002, ApJ, 580, 317 
Carter, C, Karovska, M., Jerius, D., Glotfelty, K., & Beikman, S. 2003, 

Astronomical Data Analysis Software and Systems XII; eds. H. E. Payne, 

R. I. Jedrzejewski, & R. N. Hook, 477 
Cash, W. 1979, ApJ, 228, 939 
Cohen, M., & Fuller, G. A. 1985, ApJ, 296, 620 

Curiel, S., Rodriguez, L. F, Moran, J. M., & Canto, J. 1993, ApJ, 415, 191 
Dougados, C, Cabrit, S., Lavalley, C, & Menard, F. 2000, A&A, 357, L61 
Ducourant, C, Teixeira, R., Perie, J. P., et al. 2005, A&A, 438, 769 
Dupree, A. K., Brickhouse, N. S., Smith, G. H., & Strader, J. 2005, ApJ, 625, 
L131 

Dutrey, A., Guilloteau, S., Duvert, G, et al. 1996, A&A, 309, 493 
Eisloffel, J., & Mundt, R. 1998, AJ, 115, 1554 

Favata, E, Fridlund, C. V. M., Micela, G, Sciortino, S., & Kaas, A. A. 2002, 
A&A, 386, 204 



Garcia, P. J. V, Ferreira, J., Cabrit, S., & Binette, L. 2001a, A&A, 377, 589 
Garcia, P. J. V, Cabrit, S.m Ferreira, J., & Binette, L. 2001b, A&A, 377, 609 
Glassgold, A. E., Najita, J., & Igea, J. 2004, ApJ, 615, 972 
Grosso, N., Feigelson, E. D., German, K. V., et al. 2006, A&A, 448, L29 
Giidel, M., Skinner, S. L., Briggs, K. R., et al. 2005, ApJ, 626, L53 
Giidel, M., Briggs, K. R., Arzner, K., et al. 2007a, A&A, 468, 343 
Giidel, M., Telleschi, A., Audard, M., et al. 2007b, A&A, 468, 515 
Hamann, F. 1994, ApJS, 93, 485 

Hartigan, P., Edwards, S., & Ghandour, L. 1995, ApJ, 452, 736 

Johns-Krull, C. J., & Herczeg, G. J. 2007, ApJ, 655, 345 

Hirth, G. A., Mundt, R., & Solf, J. 1997, A&AS, 126, 437 

Kitamura, Y, Kawabe, R., & Saito, M. 1996a, ApJ, 457, 277 

Kitamura, Y, Kawabe, R., & Saito, M. 1996b, ApJ, 465, L137 

Konigl, A., & Pudritz, R. E. 2000, Protostars and Planets IV (Tucson: University 

of Arizona Press), 759 
Kwan, J., Edwards, S., & Fischer, W. 2007, ApJ, 657, 897 
Lavalley, C, Cabrit, S., Dougados, C, Ferruit, P., & Bacon, R. 1997, A&A, 327, 

671 

Lavalley-Fouquet, C, Cabrit, S., & Dougados, C. 2000, A&A, 356, L41 
Leinert, C, Haas, M., Mundt, R., Richichi, A., & Zinnecker, H. 1991, A&A, 
250, 407 

Muzerolle, J., Hartmann, L., & Calvet, N., 1998, AJ, 116, 2965 

Nayakshin, S., & Kallman, T. R. 2001, ApJ, 546, 406 

Pravdo, S. H, Feigelson, E. D., Garmire, G, et al. 2001, Nature, 413, 708 

Pravdo, S. H, Tsuboi, Y, & Maeda, Y. 2004, ApJ, 605, 259 

Pravdo, S. H, & Tsuboi, Y. 2005, ApJ, 626, 272 

Preibisch, T., Kim, Y.-C, Favata, E, et al. 2005, ApJS, 160, 401 

Pyo, T-S.,Kobayashi, N., Hayashi, M., et al. 2003, ApJ, 590, 340 

Raga, A. C, Noriega-Crespo, A., & Velazquez, P. 2002, ApJ, 576, L149 

Ray, T. P., Muxlow, T. W. B., Axon, D. J., et al. 1997, Nature, 385, 415 

Reipurth, B., & Bally, J. 2001, ARA&A, 39, 403 

Rodriguez, L. F. 1995, RevMexAA, 1, 10 

Safier, P. N. 1992, ApJ, 392, 492 

Shu, F. H., Najita, J. R., Shang, H., Li, Z.-Y. 2000, Protostars and Planets IV 

(Tucson: University of Arizona Press), 789 
Solf, J., & Bohm, K. H. 1993, ApJ, 430, L31 

Stapelfeldt, K., Burrows, C. J., Krist, J. E., & the WFPC2 Science Team. 1997, 
in Herbig-Haro Flows and the Birth of Stars, IAU Symposium No. 182, Ed. 
B. Reipurth & C. Bertout (Dordrecht: Kluwer), 355 
Takami, M., Chrysostomou, A., Bailey, J., et al. 2002, ApJ, 568, L53 
Telleschi, A., Giidel, M., Briggs, K. R., Audard, M., & Scelsi, L.2007a, A&A, 
468, 443 

Telleschi, A., Giidel, M., Briggs, K. R., Audard, M., & Palla, F. 2007b, A&A, 
468, 425 

Testi, L., Bacciotti, E, Sargent, A. I., Ray, T. P., & Eisloffel, J. 2002, A&A, 349, 
L31 

Tsujimoto, M., Koyama, K., Kobayashi, N., et al. 2004, PASJ, 56, 341 
Uchida, Y, & Shibata, K. 1985, PSAJ, 37, 515 

Vuong, M. H., Montmerle, T., Grosso, N., et al. 2003, A&A, 408, 581 
White, R. J., & Ghez, A. M. 2001, ApJ, 556, 265 
White, R. J., & Hillenbrand, L. A. 2004, ApJ, 616, 998 



